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^ ABSTRACT 

HH ■ We compare the redshifts, host galaxy metallicities, and isotropic (E^^iso) and 

i-C ■ beaming-corrected (E^) gamma-ray energy release of 16 long-duration gamma- 

I ■ ray bursts (LGRBs) at z < 1. From this comparison, we find no statistically 

i-{ ■ significant correlation between host metallicity and redshift, E^^iso, or E.y. These 

^ ■ results are at odds with previous theoretical and observational predictions of 

an inverse correlation between gamma-ray energy release and host metallicity, 

On) ■ as well as the standard predictions of metallicity- driven wind effects in stellar 

evolutionary models. We consider the implications that these results have for 



> 



en . LGRB progenitor scenarios, and discuss our current understanding of the role 

(^ . that metallicity plays in the production of LGRBs. 

O 

o 

1. Introduction 

^ . Gamma-ray bursts (GRBs) are among the most energetic phenomena observed in the 

j^ . universe. The long-duration (Tgo > 2 s) subclass of these events, LGRBs, are typically 

associated with the core-collapse of massive stars (Woosley 1993). However, a detailed 
picture of the progenitors that produce LGRBs is still unclear. Potential mass stripping 
mechanisms and engine models for LGRBs include winds (e.g. Yoon et al. 2006, Woosley 
& Heger 2006), binary mergers (e.g. Podsiadlowski et al. 2004, 2010; Fryer & Heger 2005) 
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and magnetars (e.g. Wheeler et al. 2000, Burrows et al. 2007, Uzdensky & MacFadyen 2007, 
Bucciantini et al. 2009, Nagataki 2009). 

The most commonly-cited progenitor model for LGRBs was originally proposed by 
Woosley (1993), and is known as the collapsar model. Under the classical assumptions of 
this model, the progenitor is a single rapidly-rotating massive star, which maintains a high 
enough angular momentum over its lifetime to generate a LGRB from core-collapse to an 
accreting black hole. The collapsar model has been considered somewhat paradoxical since 
the progenitors, assumed to be Wolf-Rayet stars, are expected to shed their outer H and He 
shells during periods of high wind-driven mass loss. This in turn would lead to a loss of the 
angular momentum that is critical for GRB production in the collapsar model. However, 
multiple models have shown that massive stars evolving in low-metallicity environments 
retain sufficient angular momentum to produce a GRB, due to the effects of metallicity on 
line-driven mass loss rates (e.g. Hirschi et al. 2005, Yoon et al. 2006, Woosley & Heger 2006). 

In recent years, observations of LGRB host galaxies have generally supported the claim 
that the progenitors may favor low-metallicity environments. The host galaxies are generally 
fainter, more irregular, and have lower oxygen abundances than the general galaxy population 
out to 2; ~ 1 (e.g. Stanek et al. 2006; Fruchter et al. 2006; Kewley et al. 2007; Wainwright et 
al. 2007; Modjaz et al. 2008; Kocevski et al. 2009; Levesque et al. 2010a, 2010b). A number 
of studies have even proposed the idea of a "cut-ofF' maximum host metallicity for LGRB 
production (Stanek et al. 2006, Wolf & Podsiadlowski 2007, Modjaz et al. 2008, Kocevski 
et al. 2009), although the discovery of several high-metallicity LGRB host environments 
suggest that a low-metallicity cut-off is unlikely (Levesque et al. 2010b, 2010c; Graham et 
al. 2009). The high-metallicity host galaxy of the relativistic SN 2009bb presents additional 
evidence that the central-engine-driven relativistic explosions thought to power LGRBs are 
not restricted to low-metallicity progenitors (Levesque et al. 2010d, Soderberg et al. 2010a). 
Due to the young ages assumed for the massive progenitors of LGRBs (< 10 Myr; Woosley 
et al. 2002), these measured environmental metallicities can be considered representative of 
the natal abundance properties of the progenitor stars (see Simon-Diaz et al. 2006, Hunter 
et al. 2007). 

Even with these careful studies of LGRBs and their host environments, the role that 
metallicity plays in progenitor evolution and GRB production remains unclear. If metallic- 
ity does indeed have a direct impact on progenitor properties (such as angular momentum) 
that are key to producing these high-energy core-collapse events, we would therefore expect 
metallicity to have some correlation with the explosive properties of LGRBs. A low metal- 
licity environment produces stars with higher helium core masses and faster rotation rates, 
leading to LGRBs that are expected to be more energetic (e.g. MacFadyen & Woosley 1999); 
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from these assumptions low metallicities should produce LGRBs with a higher energy release 
in the gamma-ray regime {E^). 

Several previous studies have investigated this possibility. Ramirez-Ruiz et al. (2002) 
found a tentative positive correlation between the isotropic energy release {E^^iso) and the 
offset of a GRB from the center of its host galaxy (ro). This offset correlation was proposed 
to be a potential artifact of a correlation between E^^iso and low metallicity - chemical 
abundance gradients have shown that stars at higher tq, in the outskirts of their hosts, have 
lower metallicities on average (e.g. Zaritsky et al. 1994, van Zee 1998, Henry & Worthey 
1999). Since E^^iso is calculated assuming a quasi-spherical GRB explosion geometry, rather 
than accounting for the expected effects of a potential conical geometry with a narrow 
opening angle for the GRB jet {9j] see Frail et al. 2001), this result suggested that low 
metallicity was associated with either higher E^ or narrower GRB jets. 

A comparison between host metallicity and E^^iso for 5 nearby GRBs (z < 0.3) was 
performed by Stanek et al. (2006). These authors found a correlation between the two quan- 
tities. All but one of the bursts in this sample were "sub-luminous" LGRBs, a potentially 
unique class of GRB with E^^iso values that are much lower than the general population 
(e.g. Soderberg et al. 2004a, Soderberg 2006). Stanek et al. (2006) argued that this corre- 
lation supported the idea of a "threshold" metallicity for producing "cosmological" GRBs 
with more typical luminosities, given that the burst with the highest E^^iso in the sample 
was produced in the host galaxy with the lowest metallicity. However, they also cautioned 
that this conclusion was speculative due to the small size and sub-luminous nature of their 
sample. 

Wolf & Podsiadlowski (2007) also investigated the possibility of a trend relating E^^iso 
and metallicity, performing the same comparison described in Stanek et al. (2006) and in- 
cluding a sample of 13 "cosmological" {z > 0.2) LGRBs. They concluded that, while the 
Stanek et al. (2006) relation holds true for the nearby sub-luminous LGRBs, no relation 
is apparent in the larger sample. However, the metallicies derived for these "cosmological" 
LGRBs were extrapolated from the general luminosity-metallicity relation for star-forming 
galaxies, a relation that LGRB host galaxies are now known to not follow (Kewley et al. 
2007, Modjaz et al. 2008, Levesque et al. 2010a). 

To perform a robust test for a correlation between metallicity and the gamma-ray energy 
release, a large and uniform sample of LGRBs with known host metallicities, E^^iso, and 9j 
is required. Here we present the results of such a comparison, using LGRB host metallicities 
from Levesque et al. (2010a, 2010b) and E^^iso and 6j measurements for these LGRBs drawn 
from the literature. 
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2. Comparing Host Metallicity and Gamma- Ray Energy Release 

We have compiled redshifts, host metallicities, E^^iso, and 6j values for a sample of 
16 z < 1 LGRBs. The redshifts were determined from GCN circulars and other current 
literature (see Table 1). Metallicities for all but one of host galaxies were taken from Levesque 
et al. (2010a,b); for GRB 100316D we use the metallicity from Chornock et al. (2010). All 
metallicities in our sample were determined using the R23 metallicity diagnostic presented in 
Kewley & Dopita (2002) and refined in Kobulnicky & Kewley (2004) and Kewley & Ellison 
(2008). For the host galaxies of GRBs 980703 and 020405, we include both upper and lower 
metallicities since the degeneracy in the R23 diagnostic could not be resolved (see Levesque et 
al. 2010a,b). Since the majority of the hosts in our sample are unresolved, these global host 
metallicities are adopted as a proxy for progenitor metalhcity in out analyses; it is expected 
that, in general, a low-metallicity progenitor would be associated with a low-metallicity 
host. It is also worth noting that most LGRB host galaxies are low-luminosity compact 
dwarfs, which are known to have small metallicity gradients (however, see also the analysis 
of the unusual and well-resolved GRB 980425 host galaxy presented in Christensen et al. 
2008). Our metallicities also differ in many cases from those given for the nearby LGRB 
hosts presented in Stanek et al. (2006); however, these metallicities were determined from 
a variety of different spectra (Sollerman et al. 2005, Hammer et al. 2006, Gorosabel et al. 
2005, Modjaz et al. 2006) that included differing resolutions, S/N, contributions from the 
underlying supernovae, and large uncertainties in flux calibrations. We consider the results 
of our uniform host galaxy survey to be more robust. 

For all but two of the LGRBs in our sample, values for E^^iso were taken from Amati 
(2006) and Amati et al. (2008). In the case of GRB 050826, we adopt the value for E^^iso 
from Butler et al. (2007), and for GRB 100316D we adopt E^^i^o from Starling et al. (2010). 
Amati (2006) and Amati et al. (2008) determine E^^iso ranging over (10 - 10,000)/(1+^) keV, 
extrapolating from fits to data drawn from various instruments. While Amati (2006) note 
that this approach could produce systematic errors, the data are internally self-consistent. 
Butler et al. (2007) determined E^^i^o in the 1-10,000 keV band for GRB 050826. Amati 
(2006) estimates that the difference inherent in integrating from 1-10,000 keV rather than 
10-10,000 keV integrations is typically on the order of 3-5%, and no larger than 10%. Starling 
et al. (2010) find a lower limit for E^^iso in the 1-160 keV range for GRB 100316D {Swift 
satellite observations were temporarily halted before 7-ray emission from the burst had fully 
ceased). They note that there is very little flux observed for this nearby {z = 0.059) sub- 
luminous GRB above 50 keV; as a result, additional contribution to the isotropic energy 
from the higher-energy regime is expected to be negligible. 

To convert from E^^iso into E^, we adopted values for 9j taken from the literature. The 
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jet opening angle 9j is typically calculated based on the time of an observed "break" in the 
afterglow lightcurve and using the formula of Sari et al. (1999); in cases where the jet break 
time is unclear, a lower or upper limit can be placed on 9j. For GRB 030528, we adopt 
the Rau et al. (2004) approximation of a jet break at ~3 days, and use this to calculate a 
6j = 0.14 radians. For bursts that have shown evidence of a quasi-spherical, rather than 
conical, explosion geometry (GRBs 980425, 020903, 031203, 060218, and 100316D), E^^iso 
and E^ are roughly equivalent. 

We plot host metallicity against redshift, E^^iso, and E^ in Figure 1. We have included 
the redshift comparison to illustrate any potential correlation that may appear as an artifact 
of metallicity evolution with redshift. For host galaxies with both upper and lower branch 
metallicities from the Kobulnicky & Kewley (2004) R23 diagnostic, we plot two data points 
connected by a dotted line to illustrate their correspondence to a single host. We find that 
there is no statistically significant correlation between metallicity and redshift (Pearson's 
r = 0.10, p = 0.71 assuming lower-branch metallicities, and Pearson's r = 0.28, p = 0.29 
assuming upper-branch metallicities). 

Similarly, when comparing metallicity and E^^iso-, we again find no statistically significant 
correlation (Pearson's r = 0.08, p = 0.78 assuming lower-branch metallicities, and Pearson's 
r = 0.10, p = 0.72 assuming upper-branch metallicities). This result is at odds with the 
inverse correlation between E^y^iso and metallicity proposed by Stanek et al. (2006). Finally, 
we find no statistically significant correlation between metallicity and E^ (Pearson's r = 
0.18, p = 0.60 assuming lower-branch metallicities, and Pearson's r = 0.37, p = 0.26 assuming 
upper-branch metallicities). Therefore, even taking beaming effects into account we find no 
evidence for a relation between host metallicity and the true gamma-ray energy release in 
LGRBs 



3. Discussion 

The complexities of the role that metallicity plays in the production of LGRBs and the 
evolution of their progenitors can be summarized by several key findings: 

1. From a comparison of 16 2; < 1 LGRBs in this work, we find no evidence for a corre- 
lation between host metallicity and gamma-ray energy release, considering both E^^iso 
and E^. This conclusion does not agree with the predictions and claimed relations of 
past work (e.g. MacFadyen & Woosley 1999, Ramirez-Ruiz et al. 2002, Stanek et al. 
2006). This is also at odds with the standard predictions of metallicity-driven wind 
effects in LGRB progenitor evolutionary models. 
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2. There does not appear to be a clear cut-off host metalhcity for LGRBs, as has been 
previously proposed (e.g. Wolf & Podsiadlowski 2007, Modjaz et al. 2008, Kocevski et 
al. 2009); we have observed LGRBs in high-metallicity environments (Graham et al. 
2009; Levesque et al. 2010b, 2010c). There does not yet appear to be a clear maximum 
metallicity above which LGRB progenitors cannot be formed. 

3. A number of studies demonstrate that LGRBs do preferentially occur in galaxies with 
lower metallicities than the general population. LGRB hosts observed out to 2; ~ 1 
fall below the luminosity-metallicity and mass-metallicity relations for the general star- 
forming galaxy population (e.g. Stanek et al. 2006, Kewley et al. 2007, Modjaz et al. 
2008, Levesque et al. 2010a,b). However, the physical phenomena driving this trend 
remain unclear. 

The key assumptions that are refuted by items 1 and 2 - a correlation between lower 
host metallicities and higher gamma-ray energy releases, and a proposed upper metallicity 
cut-off for LGRB host galaxies - are based on the traditional collapsar model and current 
assumptions regarding the effects of metallicity on massive star evolution. Specifically, under 
the assumption of the collapsar model, lower-metallicity host environments are expected to 
produce progenitors with higher angular momentum, which should in turn produce LGRBs 
with higher gamma-ray energy releases (e.g. MacFadyen & Woosley 2001, Hirschi et al. 2005, 
Yoon et al. 2006, Woosley & Heger 2006). However, current evolutionary models for massive 
stars do not properly address the difficulties of modeling mass loss mechanisms, which may 
be anisotropic (Meynet & Maeder 2007) and could also include complex effects such as wind 
clumping (Crowther et al. 2002) and rotation-driven mass loss effects (Meynet & Maeder 
2000). According to Dessart et al. (2008), current magnetohydrodynamic simulations of 
low-metallicity massive stars actually produce core angular momenta that are too high to 
generate GRB-producing collapsars. Adopting complete and rigorous treatments of mass loss 
components and magnetic processes in massive stellar evolutionary models could potentially 
yield evolutionary pathways for collapsars that are not strongly dependent on a strict cut-off 
for the progenitor's natal metallicity. It is also possible that some of the variation seen in our 
comparison could be attributed to variations in the initial masses of the LGRB progenitor 
stars. We are not able to extrapolate any information about stellar progenitor masses from 
the data. However, it is clear that metallicity is not the sole determinant of E^^iso or E^ in 
LGRB production, and that additional parameters must be considered in future studies. 

Alternatively, it is also possible that these recent results regarding LGRBs and their host 
metallicities may be in better agreement with other alternative progenitor pathways, such 
as magnetar or binary scenarios. These models do not necessarily require a low-metallicity 
environment for the evolution and development of the critical mechanism that produces a 
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LGRB (see, for example, Fryer & Heger 2005; Podsiadlowski et al. 2010). However, it is not 
yet clear whether any of these possibilities can adequately address the phenemonon observed 
in item 3. We cannot yet provide a physical explanation for why LGRB hosts would only 
have lower metallicities relative to their mass, rather than metallicities that are uniformly 
low or fall below a particular threshold, although this may relate to the young progenitor 
ages and star formation histories of the hosts (see Berger et al. 2007). For the moment, it is 
clear that, while a lower-than- average host metallicity is a key component of LGRBs, the role 
of metallicity in progenitor evolution and LGRB production currently remains a mystery. A 
number of high-redshift host galaxies currently have absorption-line metallicities available in 
the literature (e.g. Prochaska et al. 2007, Savaglio et al. 2009, Rau et al. 2010); while the lack 
of satisfactory calibrations between emission-line and absorption-line metallicities precludes 
adding these high-redshift galaxies to the current sample, future extension of this work to 
include a greater redshift range of LGRB hosts would be valuable. Subsequent work in this 
area would also benefit greatly from an improved understanding, both observational and 
theoretical, of the various mechanisms that drive mass loss and impact late-type evolution 
in massive stars. 
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Fig. 1. — The metallicity vs. redshift (top), E^^iso (center), and E.y relations for our sample 
of LGRB host galaxies. The hosts have been separated into redshift bins by color in order 
to better illustrate redshift effects in these comparisons. Host galaxies with both lower- 
and upper-branch metallicities from the Kobulnicky & Kewley (2004) R23 diagnostic are 
indicated by lower and upper data points connected by dotted lines. Upper and lower limits 
are indicated by arrows. Host galaxies with both upper and lower limits on their E^ values 
are indicated by data points connected by solid lines. The five nearby LGRB host galaxies 
included in the original Stanek et al. (2006) relation are marked with outer circles. 
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Table 1. Properties of Nearby {z < 1) LGRBs 



GRB 


2 


log(0/H) + 12=^ 


E-,,,so (1052 erg) 


9j (radians) 


E^ (10^2 erg) 


References'' 


980425 


0.009 


o 4+0.2 


0.00010 ± 0.00002 


quasi-spherical 


0.00010 ± 0.00002 


1,2,3,4,5 


060218 


0.034 


8.21 


0.0053 ± 0.0003 


quasi-spherical 


0.0053 ± 0.0003 


6,2,7,8 


100316D 


0.059 


8.3 


> 0.0059 ± 0.0005 


quasi-spherical 


> 0.0059 ± 0.0005 


9,10,11,12 


031203 


0.105 


8.27 


0.010 ± 0.004 


quasi-spherical 


0.010 ± 0.004 


13,2,3,14 


030329 


0.168 


8.13 


1.5 ± 0.3 


0.45 


0.15 ± 0.03 


15,2,7,16,17 


020903 


0.251 


8.07 


0.0024 ± 0.0006 


quasi-spherical 


0.0024 ± 0.0006 


18,2,7,19 


050826 


0.296 


8.83 


0.03t°:°^ 


>0.2 


>0.0006 ± 0.0004 


6,20,21,22 


020819B 


0.410 


9.0 


0.68 ± 0.17 




<0.68 ± 0.17 


23,24,7 


990712 


0.434 


~8 4+°-2 


0.67 ± 0.13 


>0.71 ± 0.03 


>0.16 ± 0.04 


25,2,7,26, 27 


010921 


0.451 


8.24 


0.95 ± 0.10 


0.56 ± 0.04 


0.15 ± 0.04 


28,20,7,26, 29 


020405 


0.691 


8.33/8.59 


10 ± 1 


0.14 ± 0.02 


0.10 ± 0.04 


30,2,7,26, 31 


991208 


0.706 


8.02 


22.3 ± 1.8 


<0.15 ± 0.01 


<0.25 ± 0.06 


32,20,7,26, 33 


030528 


0.782 


o 4+0.2 


2.5 ± 0.3 


0.14 


0.024 ± 0.003 


34,20,7,35 


051022 


0.807 


8.62 


54 ± 5 


0.08 


0.17 ± 0.02 


36,2,7,37 


050824 


0.828 


~8.4 ± 0.2 


0.130 ± 0.029 


0.08 


0.0004 ± 0.0001 


6,20,3,38 


980703 


0.966 


8.31/8.65 


7.2 ± 0.7 


0.20 ± 0.02 


0.14 ± 0.01 


39,20,7,26, 40 



''Host metallicity, derived from the Kobulnicky & Kewley (2004) R23 diagnostic. Metallicitics have a systematic 
error of ±0.1 dex except where noted. 

''References: (1) Tinney et al. (1998), (2) Levesque et al. (2010a), (3) Amati (2006), (4) Kulkarni et al. (1998), (5) 
Li & Chevalier (1999), (6) Swift data archive, (7) Amati et al. (2008), (8) Soderberg et al. (2006), (9) Sakamoto et al. 
(2010), (10) Chornock et al. (2010), (11) Starling et al. (2010), (12) Soderberg et al. (2010b), (13) Mereghetti & Gotz 
(2003), (14) Soderberg et al. (2004b), (15) Greiner et al. (2003), (16) Berger et al. (2003), (17) Frail et al. (2005), 
(18) Ricker et al. (2002), (19) Soderberg et al. (2004a), (20) Levesque et al. (2010b), (21) Butler et al. (2007), (22) 
Mirabal et al. (2007), (23) Hurley et al. (2002), (24) Levesque et al. (2010c), (25) Galama et al. (1999), (26) Bloom 
et al. (2003), (27) Fruchter et al. (2000), (28) Djorgovski et al. (2001), (29) Price et al. (2002), (30) Masetti et al. 
(2002), (31) Price ct al. (2003), (32) Hurley & Cline (1999), (33) Sagar ct al. (2000), (34) Atteia et al. (2003), (35) 
Rau et al. (2004) (36) Hurley et al. (2005), (37) Nakagawa ct al. (2006), (38) Racusin et al. (2009), (39) Djorgovski 
et al. (1998), (40) Frail ct al. (2003) 



